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Introduction 



The notion of quasi-free state was developed in the framework of the C*-algebraic 
approach to the canonical commutation relation (CCR) [Ml [IHl E, |16]. The CCR- 
algebra is generated by the Weyl unitaries (satisfying a commutation relation, therefore 
Weyl algebra is an alternative terminology). The quasi-free states on CCR can be 
regarded as analogues of Gaussian distributions in classical probability: The n-point 
functions can be computed from the 2-point functions and in a kind of central limit 
theorem the limiting state is quasi-free and it maximizes the von Neumann entropy 
when the 2-point function is fixed [17]. The quasi-free states are quite tractable, for 
example the von Neumann entropy has an explicit expression [SI Ej . 

The Markov property was invented by Accardi in the non-commutative (or quantum 
probabilistic) setting [H O [3l Hj. This Markov property is based on a completely 
positive, identity preserving map, so-called quasi-conditional expectation and it was 
formulated in the tensor product of matrix algebras. A state of a tensor product system 
is Markovian if and only if the von Neumann entropy increase is constant. This property 
and a possible definition of the Markov condition was suggested in [TH]. A remarkable 
property of the von Neumann entropy is the strong subadditivity [ISl El UHl |19] which 
plays an important role in the investigations of quantum system's correlations. The 
above mentioned constant increase of the von Neumann entropy is the same as the 
equality for the strong subadditivity of von Neumann entropy. 

A CCR (or Weyl) algebra is parametrized by a Hilbert space, so we use the notation 
CCR(?i) when 7i is the Hilbert space. Assume that ipus is a state on the composite 
system CCR(7Yi) ® CCR(?i2) ® CCR(?^3). Denote by V5i2,v^23 the restrictions to the 
first two and to the second and third factors, similarly (p2 is the restriction to the 
second factor. The Markov property is defined as 

where S denotes the von Neumann entropy. When 99123 is quasi-free, it is given by 
a positive operator (corresponding to the 2-point function) and the main goal of the 
present paper is to describe the Markov property in terms of this operator. The paper 
[20] studies a similar question for the CAR algebra and [5] is about multivariate Gaus- 
sian distributions. Although the multivariate Gaussian case (in classical probability) is 
rather different from the present non-commutative setting, we use the same block ma- 
trix formalism (and the paper [S] was actually a preparation of this problem). The proof 
of the main result uses the description of sufficient statistics in the non-commutative 
case. A quasi-free state is described by a block matrix and the Markov property is 
formulated by the entries. A Markovian quasi-free state induces multivariate Gaussian 
restrictions, but they are very special in that framework. 

The paper is organized as follows. The preliminary section contains some crucial 
properties of the Weyl unitaries, the Fock space, the CCR algebra and quasi-free 
states. This is written for the sake of completeness, the results are known but not 
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well-accessible in the literature, cf . [6l [3, [TO] . The main point is the von Neumann 
entropy formula which is well-known for the CCR quasi-free state. In the next section 
we investigate the quasi-free Markov triplets. We obtain a necessary and sufficient 
condition described in the block matrix approach: The block matrix should be block 
diagonal. There are nontrivial Markovian quasi-free states which are not a product in 
the time localization. The existence of such state is interesting, because it is in contrast 
to the CAR case ^U\. However, the first and the third subalgebras are always inde- 
pendent. Finally we prove that commuting field operators form a classical Gaussian 
Markov triplet. 



1 CCR algebras and quasi-free states 
1.1 Introduction to Weyl unitaries 

In this part the basis of Hermite functions of the Hilbert space is described in 

details, the creation, annihilation operators and the Weyl unitaries are constructed. 

The Hermite polynomials 

Jn 

H^ix) := (-l)"e-^_e--' (n = 0, 1, . . .) (1) 

are orthogonal in the Hilbert space L^(M, e~^^ dx), they satisfy the recursion 

- 2xHn{x) + 2nHn-i{x) = (2) 

and the differential equation 

H'^ix) - 2xH'^{x) + 2nHn{x) = 0. (3) 
The normalized Hermite polynomials 

Hnix) = -^J==H^{x) (4) 



form an orthonormal basis. From this basis of L^(]R, e ^ dx), we can get easily a basis 
in L2(M): 

^n{x) := e-'^'/^Hnix) . (5) 

These are called Hermite functions. In terms of the Hermite functions equation ([2]) 
becomes 

If the operators a and a+ are defined as 



ai^n = Vn(pn-i, a+V9„ = Vn + l^Pn+i (7) 
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with ai^Q = and the multiphcation by the variable x is denoted by Q, then ([6]) is 

Q = -^(a + a+). (8) 

From the equation 
one can obtain 

P^n := = -j^ , (9) 

that is 

P=^{a-'-a). (10) 

Therefore, 

a = i=(g + iP), a+ = -L(Q_iP). 
Lemma 1.1 For z E C the identity 

e(^) := —F^fn{x) = TT^^/^exp ( — ) exp{zxV2) 

n=0 * 

/io/(is. Moreover, 

e(z) = e^'^Vo, ||e(z)|| =el^l'/l 
Proof. The identity can be deduced from the generator function 

=exp(2xt-t=^) (11) 

n=0 

of the Hermite polynomials. □ 
The above e{z) is called exponential vector. 

For 2; G C, the operator i{za — za^) is defined originally on the linear combinations 
of the basis vectors (pn and it is a symmetric operator. It can be proven that its closure 
is self-adjoint, therefore exp(2;a — za^) becomes a unitary. 

W{z) := e'"-''"* (12) 

is called Weyl unitary. Note that 

W{z) =expiV2(aP + /?g) 

ii z = a + i(3. Multiple use of the identity 

gi(tQ+«p) _ exp(itM/2)e"^e'"-^ = exp(-itM/2)e'"-^e'*^ (13) 

gives the following result. 
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Theorem 1.1 

W{z)W{z') = W{z + z')exp{ilm{zz')) 

for z, z' G C. 

With straightforward computation one gets the foUowing. 
Lemma 1.2 

Mz)\\ 

The functions 

Lni^) = E M^^^^t^^M, ^' (" > (14) 

^ - + fc)! ^ ^ ^ ^ 

are caUed associated Laguerre polynomials. We write simply for a = 0. 

Theorem 1.2 For n > m 



holds. 



Proof. First note that definition implies 

a^(a+)>o = <! (15) 
[ if k> n. 

If [A, S] commutes with A and S, then the formula e^e^ = ^{^M/'^^a+b ]^q\(\^ Since 
[—za^, za] = \z\'^I, we can write 

.12, 



-|2|2/2„-2a+ °° k 
V "'• I,— n 



fe=0 



Now we can compute the matrix elements: 

(n- A;)! 
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g " (—zYz''n\m\ 



u|2/2 (-2;)^^'' I mini 

2^ 2^ -&kry {n-ky.im-iy}'^-'' "^-'^ 



(—zYz'' I mini 



EE 

il kl \ (n — k)l (m — i)l 



Jm—Ln—ki 



where 6k/ denotes the Kronecker symbol. For n > m, we get non-vanishing elements 
if and only if k = n — m + i, where n — m < k < n and by the formula ( |T4l) we obtain 

V n! ^!(m — /)!(n — m + /)! 



n\ 



as we stated. □ 
Note that the case m > n can be read out from the Theorem, since 



{(Pm,W{z)(pn) = {(Pn,W{-z)(p^). 

The case m = n involves the Laguerre polynomials. The analogue of ffTTj) is the formula 

Et"^n(^) = Y^exp(-^) (16) 

n=0 

which holds for |t| < 1 and x G M^. This formula is used to obtain 

00 I |2 1 I 

- ^)(^„, W{z)^n) = exp ( - ^Y^) (17) 

n=0 ^ 

for < < 1. Note that 



Z^ = ^/i"(l-/i)|y.„)(y.„| (18) 

n=0 

is a statistical operator (in spectral decomposition). In the corresponding state the 
self-adjoint operator 

= \—\z)a + (— iz)a 

i 

has Gaussian distribution. 
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1.2 The Fock space 

Let 7i be a Hilbert space. If vr is a permutation of the numbers {1,2, ... ,n}, then on 
the n-fold tensor product 7i®" := H ® H <^ . . . <^ TC we have a unitary [/^ such that 

UMl ® /2 ® • • • ® /n) = fnil) ® /7r(2) ® • • • ® /ttH- 

The operator 

PnU'l ® /2 ® • • • ® /n) := ^ 5Z •^^(1) ® •^'^(2) ® • • • ® /^(„) 

is a projection onto the symmetric subspace 

7^^" := G 7^^" : f/^^ = g for every vr} 
which is the hnear span of the vectors 

l/l, /2, • • • , /n) = /l V /2 V . . . V := ^ V ® /^(2) ® . . . ® , 

where /i, /2, • • • , /n e Obviously, 

/i V /2 V . . . V /„ = /^(i) V /^(2) V ... V 

for any permutation vr. 

Assume that ei, 62, . . . , is a basis in 7i. When we consider a vector 

ei(i) V ei(2) V ... V 

in the symmetric tensor power 7i^", we may assume that 1 < < i{2) < . . . < 
i{n) < m. A vector et may appear several times, assume that its multiplicity is rt, that 
is, Tt := {i : = t}. The norm of the vector is \/ri\r2\'77~r\J and 

\ , , ] ^e,(i) V ei(2) V ... V ei(„) : 1 < < t{2) <...< i{n) < m} (19) 

is an orthonormal basis in T-^^". Another notation is 

|e?, 62', . . . , e"^) = ei(i) V 6^(2) V ... V ei(„). 

The symmetric Fock space is the direct sum 

where $ is called the vacuum vector and in this spirit the summand Ti"^" is called 
the n-particle subspace. Since 7i^^ is identical with 7i, the Hilbert space jF(7i) is an 
extension of Ti. The union of the vectors ( fT9l) (for every n) is a basis of the Fock space. 
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Lemma 1.3 IfU^Hi® Hi, then J'{H) = T{Hi) ® T{H2). 
Proof. It is enough to see that 

If ei, 62, . . . , is a basis in TYi and /i, f2, ■ ■ ■ , fk is a basis in 7^2, then the (non- 
normalized) basis vector 

ej(i) V ei(2) V ... V ei(t) V V /j(2) V ... V fj{n-t) 

can be identified with 

ej(i) V ej(2) V ... V ej(t) O /,(!) V /j(2) V ... V fj(n-t) 

which is a basis vector in 7i^* . □ 

For f & H the creation operator a"'"(/) is defined as 

a+(/) l/l, /2, . . . , /n) = I/, /l, /2, . . . , /n). (20) 

a'^if) is linear in the variable / and it maps the n-particle subspace into the (n + 1)- 
particle subspace. Its adjoint is the annihilation operator which acts as 

n 

«(/) l/l, /2, ■ ■ ■ , fn) = fi)\fl, fi-1, fi+1, ...Jn). (21) 

i=l 

Given an operator A G B{7i) acting on the one-particle space, we can extend it to 
the Fock space as follows. 

n 

J^(A)\f^, /2, ...,/„) = J] l/i, . . . , fi-1, Af,, ...,/„). (22) 

1=1 

The next lemma can be shown by simple computation. 
Lemma 1.4 For f,g eH, we have 

J^{\f){g\)^a^{f)a{g). 

Another possibility for extension, or second quantization, of an operator U e B{H) 
is given by 

T{U)\f^, /2, . . . , fn) = \Ufi, Uf2, . . . , (23) 

It is easy to see that 
Lemma 1.5 

T{U^U2) = T{U^)T{U2) and T{U*) = T{Uy. 
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If ?7 is a unitary, then T{U) is unitary as well. Moreover, if U{t) is a continuous one- 
parameter group with generator A, then the generator of the continuous one-parameter 
group T{U{t)) on jF(7i) is the closure of J^{A). To show an example, we note that the 
statistical operator ( fTSjl is (1 — /i)r(/i) in the case of a one- dimensional TC. 

Lemma 1.6 Let 71 = Hi® 712 and U = Ui®U2- Then 

T{Ui®U2) = T{Ui)®T{U2). 

1.3 The algebra of the canonical commutation relation 

Let 7i be a finite-dimensional Hilbert space. Assume that for every / G 7i a unitary 
operator W{f) is given such that the relations 

WUi)WU2) = WUi + /2) exp(i /s)) (24) 

W{-f) = W{fr (25) 

hold for /i, /2, / G with (t(/i, /2) := Im(/i, /2). The abstract C*-algebra generated 
by these unitaries is unique and denoted by CCR(?-^). The relation shows that 
W{fi) and W{f2) commute if /i and /2 are orthogonal. Therefore for an n-dimensional 
H, the algebra CCR(7Y) is an n-fold tensor product 

CCR(C)®...®CCR(C). 

Since W{tf)W{sf) = W{{t + s)f) for t, s G M, the mapping t ^^ W{tf) is a one- 
parameter unitary group which is not norm continuous since ||W^(/i) — Vr(/2)|| > 
when /i ^ /2 [16]. 

The C*-algebra CCR(7i) has a very natural state 

a;(Ty(/)):=exp(- 11/1172) (26) 

which is called Fock state. The GNS-representation of CCR(7i) is called Fock rep- 
resentation and it leads to the the Fock space J-'{T-C) with cyclic vector If /i and 
/2 are orthogonal vectors, then 

uiWih)Wif2)) = a;(iy(/i + /2)) = exp(-||/i + /2||V2) 

= exp ( - II AIV2) exp ( - ||/2f /2) = u;{W{h))uj{W{f2)). 

Therefore u is a product state and it follows that the GNS Hilbert space is a tensor 
product. (This is another argument to justify Lemma IL3[ ) We shall identify the 
abstract unitary W{f) with the representing unitary acting on the tensor product 
GNS-space TiTi). The map 
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is an so-continuous 1-parameter group of unitaries, and according to the Stone theorem 

MW{tf))=expiitB{f)) 
for a self-adjoint operator B{f), called field operator. Let 



Then 



B^{f) = -{B{f)Tmf))- 



[B-{f),B^g)] = {g,f) 



is the canonical commutation relation for the creation operator B^{g) and the an- 
nihilation operator B~{f). When H = then 



W{z) = expi(a(z) + a^{z)), 



where a'^(z) = iza'^. 



1.4 Quasi-free states 

The Fock state (l26l) can be generalized by choosing a positive operator A G BiTi): 

u:A{W{f)) := exp ( - ||/f /2 - (/, Af)). (27) 
This is called quasi-free state [Ij]. 

Example 1.1 Assume that Ti = Tii © 7^2 and write the positive mapping A G B{T-C) 
in the form of block matrix: 

_A2l A22 

If / G then 



A 

1^21 A2 

ooAiWif © 0)) = exp ( - 11/11 V2 - (/, Anf)). 
Therefore the restriction of the quasi-free state ua to CCR(7ii) is the quasi- free state 

If 

^^Un 

[ A22J ' 

then Ua = ^Ah ® ^^22 • ^ 



Example 1.2 Assume that 7i is one-dimensional and let A = A > 0. We can read out 
from formulas (fTT]) and (fT^ that the statistical operator of ux in the Fock representa- 
tion is 



n=0 



+ A 
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(Note the = A/(l + A) in ([T8]).) Moreover, 

uxia+a) = A. (29) 

One can easily compute the von Neumann entropy of the state oox from the eigenvalues 
of the statistical operator Dx: 

^(a;,)= 77(A) -r/(A + l), (30) 

where 77(A) = —A log A. 

The case of finite-dimensional Ti. can be reduced to the one-dimensional by the 
spectral decomposition of the operator A. □ 

Assume that is a state of CCR(7i). If 

can be defined, then it will be called 2-point function of u. 

Theorem 1.3 Assume that the spectral decomposition of < A E B{7i) is 

m 

A = J]A,|ei)(e,|. (31) 

i=l 

Then the statistical operator of the quasi-free state uja in the Fock representation is 



^TqV)E^^(TT^)'■';^)|^-^?•■■■-">(e^e^...^ (32) 



,1=1 / r, \i=l 



where summation is over n = 0,1,2 .. . and the decompositions n = ri + r2 + . . . + Vm- 
Moreover, 

UAia^if)aig)) = {g,Af) if,gen) 

and 

S{uja) = Tt7]{A) ~Tt7]{A + I). 

Proof. The basic idea is the decomposition 

= t^Ai ® t^A2 ® • • • ® ^^\m (33) 

when the space Ti. is decomposed into the direct sum of the one-dimensional subspaces 
C|ej) and J-'{T-C) and CCR(7i) become tensor product. The statistical operator of ujx. 
is 

-^.-Et^(t^)";^K')«i 

ri=0 

the tensor product is exactly the stated matrix. 
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When we want to check the 2-point function, it is enough to consider the case 
f = g = Ci. This is OK due to 



The von Neumann entropy is deduced from ( !30l) and ( !33ll . □ 
If (ED holds, then 



■^1 5 "^2 5 • • • 5 "^m 



r(A(J + A)-')\e[\el\ . . . , e-) = n ( 

i=l ^ 

and we have 

Da = —T(A(I + A)-^) , where ca = Tr r(A(/ + A)"^). (34) 
This leads to the following result. 

Theorem 1.4 Let uja and ujb be quasi-free state o/ CCR(7i) which correspond to the 
operators < A, B & B{T-l). Their Connes cocycle is 

[DuA, DuB]t = UtT(^{A{I + A)-'f{B{I + ByY'') (35) 



where 



ut= (^Trr(A(J + A)^^)) '*(^Trr(B(j + s)-i; 



Theorem 1.5 Letu be a state ofCCR{TC) such that its 2-point function is Lj{a^{f)a{g)) 
{g, Af) (f,g&TC) for a positive operator A e B{T-t). Then S{lj) < S{uja) and equality 
implies u = uja- 

Proof. Consider the one- dimensional case when A = X. We compute the relative 
entropy S{uj\\ujx)- 



S{u!\\uJx) = —S{uj)—uj(\ogD 



A 



= -S{UJ) -log(l + X)uj(^^\ipn){<fn\^ - \0gY^^(^^n\ipn){^n 

n=0 n=0 

= -S{u;) - log(l + A) - A log = -S{ij) + SM 

i + A 

Since the relative entropy S{uj\\u!x) > if u and ux are different, the statement is 
obtained. 

The general case can be proved by similar computation. The result was also obtained 
in connection with the central limit theorem [TTI. □ 



oo 
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2 Markov triplets 



Let Ti. = Til © 7^2 © '^a be a finite-dimensional Hilbert space and consider the Fock 
representation of CCR(7^) = CCR(Hi © 7^2 © ^s) on © 7^2 © ^s)- Instead of 

the C*-algebra, we work with the weak closure in the Fock representation: ^123 := 
B{T{rii © © ^3)) = BiJ^iHi)) © B{J^{n2)) © BiJ^iHs))). This algebra has sub- 
algebras 

Au- = B{T{ni®n2®Q))~B{T{ni))®B{T{n2))®ci, 

A3: = B{J^{Q®'H2®'H^))=CI®B{J^{1-L2))®B{J^{T-Lz)), 

A2: = B{j^{o®n2®o)) = ci®B{j^{n2))®ci. 

Assume that D123 is a statistical operator in ^123 and we denote by D12, -D25 -D23 
its reductions in the subalgebras ^12,-4.2,^3, respectively. These subalgebras form a 
Markov triplet with respect to the state D123 if 

S{Du3) - S{D2s) = S{Du) - S{D2), (36) 

where 5* denotes the von Neumann entropy and we assume that both sides are finite 
in the equation. The state uj corresponding to the statistical operator -D123 is called 
Markov state. 

Condition fl36l) is equivalent to several other conditions, see, for example. Chapter 9 
of [19\ about the details and proofs. In most studies about the strong subadditivity of 
the von Neumann entropy and the equality case (l36ll . the Hilbert space is assumed to 
be finite-dimensional. In our setting the Fock space is always infinite-dimensional, so 
[12] might be the optimal reference. Here we prefer equivalent formulation in terms of 
Connes cocyles: 

[^^^123, D{lpi ® uji2)]t = [D{uJi2 ® ^3), D{ipi ®uj2® ^53)]* (37) 
for every real t, where ipi and ip^, are arbitrary states. 
Let 



An 


A12 


Ai3 


A21 


A22 


A23 


A31 


A32 


A33 



be a positive operator on Tii ®H2®'Hs and assume that D123 is the statistical operator 
of the quasi- free state ua- oja restricted to ^23 is a quasi-free state induced by the 
operator 

A22 A23 
_ A32 A33 _ 

Set 





'I 










'Au 


A12 


0" 




'I 





0" 


D = 





A22 


A23 


B = 


A21 


A22 





, c = 





A22 










A32 


A33 










/ 










/ 
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Then equality fl37jl may become 

[Dua, DujD\t = [Dujb, D^^c]t- 
According to Theorem 11.41 this is the condition 



(39) 



with a set of numbers Xt- 

One can see from formula ( |23l) that T{U) = XT{V) for a A G C implies A = 1 and 
U = V. Therefore condition (1371) becomes the following. 



Theorem 2.1 For a quasi-free state ua the Markov property ^3m) is equivalent to the 
condition 

A^\I + A)-'''D-''\I + Df = B'\I + B)-'''C-'\I + Cf (40) 
for every real t. 



The problem is the solution of this equation. Note that if condition 
every real t, then analytic continuation gives all complex t. 



holds for 



Corollary 2.1 If A gives a Markov triplet, then U*AU gives a Markov triplet as well 
when U = Diag{Ui, U2, U3) with unitaries Ui G B(Hi), 1 < i < 3. 



Example 2.1 The following matrix satisfies condition (140 



A 





All 


[a 


0] 











'a*' 




c 


0" 




'0" 














d_ 




_b 












b* 


A33 





All ci 

a* c 





d b 
b* As3 



(41) 



where the parameters a, b, c, d (and 0) are matrices. This is a block diagonal matrix, 
A = Diag(Ai,A2), so we have f{A) = Diag {f {Ai) , f [A2)) for a function /. The 
matrices B, C and D are block diagonal as well: 

B = Diag(A,Diag(rf,/)), 

C = Diag (Diag (/, c) , Diag (d, /) ) , 

D = Diag(Diag(J,c),A2). 



Therefore, 
and 



f{A)g{D) = Diag(/(Ai)Diag ((?(/), ^7(c)),/(A2)^7(A2)) 



fiB)giC) = Diag(/(Ai)Diag ((?(/), ^?(c)),Diag(/(rf)(7(rf),/(/)(7(/))). 
If fg = l,thenf{A)g{D) = f{B)g{C). 
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Note that 



An 











A22 


^23 





^32 


^33 



An 


A12 





A21 


A22 











^33 



A= Aoo Ao-^ and A 



are particular cases. 

On the basis of the previous Corollary we can use block diagonal unitaries to have 
further examples. □ 



Theorem 2.2 The condition 

A-\I + A)D{I + D)-^ 

implies that 



A 
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B-^{I + B)C{I + C)-^ 
AuA^2^A23 and A13 = ^12(^22 + iy^A23. 



(42) 
(43) 



Proof. In the computation of the inverse of a block matrix, the following formula is 
very useful. If P and S are square matrices and S is invertible, then 



1 _ 


' P Q' 


-1 


{M/sy^ 


-{M/sy^Qs-^ 




R S 






f s-^R{M/sy^Qs-^ 



(44) 



where (M/S) = P — QS ^R is the Schur complement of S* in M (see [H], Section 7.7], 
actually, the checking is a simple multiplication). If P is also invertible, the equation 



s-^ + s~^RiM/sy^Qs-^ = (M/py^ 



(45) 



also holds, where (M/P) = S — RP is the Schur complement of P in M. For 
solving fl42p . we partition the block matrix A in the following way 





'An 


A12 


A 


A = 


A21 


A22 


A 




Au 


A32 


A 



P Q 

Q* S 



where we used the fact that A is positive self-adjoint and used the notations 

^22 ^23 



S 



A*23 A33 



P = An and Q = [Au A13]. With the help of we get 

A~\I + A) = 
Similarly we write D in 2 x 2 matrix form: 



i + {A/sy' 
-s-^Q*{A/sy^ 



(A/sy^Qs-^ 
i + iA/py^ 



(46) 





'I 








D = 





^22 


^23 







^32 


^33 



I 

s 
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where denotes the block matrix [0 0], or its transpose. Now we can compute the 
left-hand side of (02]): 



A-\I + A)D{I + D)-^ = 
With a similar procedure we have 

B 



^[I + iA/S)-'] 
'.S-'Q*{A/S)-' 



and 



C 



An 


A12 





A21 


A22 











I 


"/ 











A22 











/ 



-{A/S)-^Q{I + S)-^ 
[I + {A/P)-^]S{I + S)-^ 



P Q 

Q* S 



I 

s 



(47) 



where 



S 



A22 
I 



Q = [Ai2 0] 



and the same remark concern for the block matrices as above. We get for the right- 
hand side of (1121): 



B-\I + B)CiI + C) 



-1 



I + iB/S)-' 
\S-^Q*{B/S) 



-(B/S)-'Q{I + S)-' 

[i + {B/p)-^]s{i + sy 



■ (48) 



From the equality between (l47j) and (148|) we have equations for the block matrices. 
The equality of (1,1) elements implies (A/S) = (B/S). This and the equality of 
(1,2) elements gives the equation Q{I + S)^^ = Q{I + S)^^, which lead us to A13 = 
Ai2{A22 + I)~^A23. From the (2, 1) elements we have S^'^Q* = S^^Q*, this implies the 
other necessary condition. The (2, 2) elements will be equal automatically when these 
conditions hold. □ 



According to [5j 



means that the (1,3) element of A ^ and {A + I) ^ are 0. It 



is interesting that if we take the determinant of equation 



then we have 



(Det A) (Det C) (Det D)-\DetB)-^ 

= (Det (/ + A)) (Det (/ + C))(Det (/ + D))-^(Det {I + B))-\ 

According to Theorem 5 in [5], both sides are smaller or equal than 1 and (H3l) is 
equivalent to the condition that both sides are exactly 1. 

Let X be the inverse of the block matrix 
computation yields that 



and suppose that (H2i) holds. Tedious 



Xn = {An- AuA2^A2i) \ 

X12 = — {An — A12A22 A21) A12A22, 
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Xis = 0, 

X22 = (^22-^21^11^12) ^ + ^22^^23(^33-^32^22^23) ^A32A^2 

X23 = —A22A23 (A33 — ^32^22^^23) , 

X33 = {A33- A32A^^A23y' . 



The next example shows that conditions (H3l) are not sufficient, in contrast to the 
classical Gaussian Markov triplets [5]. 



Example 2.2 The matrix 



A 



4 
5 



1 1 

_2 _2 
7 7 



14 
_ J_ 

49 



14 
_ J_ 

49 



"1 -l" 




'6 


0" 




" 1 1 " 


.1 -1. 







3 




2 2 
. 7 7 . 



1_ 

L4 
]_ 

L 14 



14 
1 



_L 

4C 
J_ 

49 J 



49 
1 



2 

■7 
2 

' 7 J 



is positive and fulfills conditions ( H3l) . but ( HOI) does not hold. Indeed, numerical 
computation shows that 

log A(/ + A)-' + logC(/ + C)-^ ^ log B{I + BY^ + log + D)-\ 

or an alternative argument is that the matrix is different from (HTl) . cf. Theorem 12.31 
This example shows that condition 

^123^2"3' = ^12^2"' (49) 

is weaker than 



-^^123 -^23 



(50) 



Note that in the finite-dimensional case ( l50l) is equivalent to 



^1/2^-1/2 

-^123 -^^23 



^1/2^-1/2 
-^^12 ^2 ' 



see Chapter 9 of 

In the notation 

K := A-\I + A), L := D'^I 
condition (HOl) becomes 



□ 



D), M := B-\I + B), N := C-\I + C) 



Theorem 2.3 The Markov property 
P G B{l-L) such that P\ni = I, P\H3 
diagonal in the form 



(51) 



is satisfied if and only if there is a projection 
and PA = AP. In other words, A is block 
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Proof. We write the matrices of the relation fl5T]) in block form: 





' Ku 


K12 


' 




' 2 








K = 




K22 


K23 


, L = 





L22 


L23 






K32 


K33 







L32 


L33 




' Mn 


M12 


" 






■ 2 








M = 


M21 


M22 







N = 





N22 













2 












2 



Suppose that the Markov property is satisfied, and we use it in the form (!5T]) . Since 
K = I + A~^, the block diagonal structure fl4ip of A is equivalent to that property of 
K. We shall work on K. 

Let C be the subalgebra generated by the set {i^^^L^'* : t G M}. By the factorization 
result in [12], there are positive matrices X,Y E C and < Z E B{T-Ci ®'H2® 'H3), 
such that 

K = XZ, L = YZ, ZX = XZ, ZY = YZ. (52) 
Since (EI]) implies that C C B(Hi © 0.2) © CI, we have X and Y in the above form 



X :-- 



X 
/ 



We write Z in a similar block form 

Z 



and Y 



Z z 
z* Z33 



Y 
/ 



where Z G 5(7^i © 7^2) and z* = [Z31, Z32]. Then 



K = XZ 
This implies that 



XZ Xz 

z* Z33 



and L = YZ 



YZ Yz 

z* Z33 





L23 



Xz = Yz = z-- 
Z commutes with X and Y and 

[K3i,K32,K33] = [Z*,Z33] = [0,^32,^33] 

In particular, K31 = K13 = and K23 = L23. 
By ([53]) and (154]), we have 



Ku Ku 
K21 K22 





L23 



XZz = ZYz 



' 2 













L22 


z = 


L22L23 



(53) 



(54) 



and we get K12K23 = K12L23 = 0. 
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If the range of L23 is 7-^2, then K12 
cases K is block diagonal. 



0, and if L 



23 



0, then 



23 



0, SO in both 



Suppose now that the range of L23 is not 7^2 and L23 7^ 0. Then there is a decom- 
position H2 = JCa ® JCb, where )Cb is the range of L23. Next we work in the frame of 
the decomposition {Hi © K-a) © K,b- 



For each vector ^ G /C 



we have = = ^. It follows that there are matrices 



JCa), such that 



X 



' Xi 


" 




' n 


" 





/ 


, Y = 





/ 



If we write 



then 



Zi 

4 



xz 



Z33 



XiZi 

4 



Zl e B{Hi © Ka) and Z33 G B{lCb 



Xizi 
Z33 



Kii 


KI2 




Kli 


KI2 


^23 


Kli 


KI2 


^3^3 



and 



YZ 



YiZi 



Yizi 
Z33 



2 











ri 
-^22 


ri 

-^^23 





ri 

-^^32 


ri 

-^^33 



with the block decompositions of the matrices in B(Hi © /Ca © 
sides. This implies that Zi commutes with both Xi and Yi, K33 = 





on the right-hand- 

■^33 



L33 and 



XiZi — YiZi 



Zl 



ri 

-^23 



(55) 



and we are in a similar situation as before. (Compare with the relations (152!) . (l53l) and 
(13^ .) We also get iiT^g-^L — exactly as before. Note that now we can write 





' Kii 


KI2 








K = 


Kli 



KI2 
ri 

-^32 


ri 

-'^23 

ri 

-'^33 




-^23 










L32 


-^33 



Again, if L23 = or if the range of L23 is /Cq, then the matrix K is block diagonal. If 
this condition does not hold, then the above procedure can be repeated: we decompose 
the subspace /Cq = /C^ © /C^, where /C^ is the range of L23, and write Xi, Yi in the 
block- diagonal form, using ( l55l) . etc. 

After repeating this procedure ra-times, we get the matrix K in the form 



K 



Kii 


-"■12 













-"■22 


rn 
-^23 











TV. 

^32 


rn 
^33 


T n' 
-^34 











T n' 
^i3 


T n' 
-^44 


-^23 











L32 


-^33 
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in B{Hi®lC^®lC'i®lC^' ®H^), where /C^' = ®ll\lCl, fCl = JC^ and i^i^a^^a = 0- Since 
7^2 is finite-dimensional, there must be some n, such that the matrix is either or 
has range /C^. In both cases, the matrix K has a block diagonal form, and so does the 
matrix A. □ 

The CCR Markov triplets have some similarity to Markov states on a product algebra 
M„(C) O (M„(C) O Mt{C)) (g) Mn{C) {n = ut). If lui is a state on M„(C) M„(C) 
and UJ2 is a state on Mt{C) ® M„(C), then ui ® u;2 is Markovian, but there are other 
Markov states, however they are constructed essentially by this idea [S]. 

3 Connection to classical Gaussians 

Markov triplets in the classical Gaussian case were studied in [S]. The present non- 
commutative situation has some relation to the classical Gaussian. 

Lemma 3.1 Let 61,62, ■■■,6^ be linearly independent unit vectors in H such that 
{6i,6j) is real, I < i < j < k. Then the Weyl unitaries W{t6j) = exp(tiB{6j)) 
commute. With respect to a quasi-free state ( [F^ , the joint distribution of the field 
operators i?(ei), -8(62), • • • , -B(efc) is Gaussian. 

Assume that /i, /2, ■ ■ ■ , /fc are orthonormal vectors and Sfi = 6i for a linear map- 
ping S. The covariance is the matrix of the linear operator S*{I + 274)5* in the basis 

fli f2i---ifk- 

Proof. The characteristic function of the joint distribution is 

(ti,t2, . . ■ ,tj) ^ uJAiexp{itiBi)exp{it2B2) . . .exp{itkBk) 

= UJAiW{ti6i+t262 + ■■■+tk6k)) 

= exp ( - UE^,,t^tM, (/ + 2A)e,))) 

= exp ( - UE^.tMSfu (/ + 2A)Sf)))). 

This gives the result. □ 

Next we assume that TC = Hi (B Ti.2 ® and assume that dimTij = (1 < i < 3). 
Choose pairwise orthogonal unit vectors fj, 1 < j < such that 

f(i-i)k+r e^i, l<i<3, 0<r<A;-l 

and unit vectors 6j, 1 < j < 3k such that 

e(i_i)fc+,, G l<i<3, 0<r<k-l (56) 

and 

(et, 6u) is real for any I <t,u < 3k. (57) 
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There is an invertible block diagonal matrix S = Diag(5'i, 5*2, 5*3) such that Sfj = ej, 
l<j< 3k. 

The Weyl unitaries W{tej) = exp{tiB{ej)) commute. The joint distribution of the 
field operators B{ej) is Gaussian with covariance block matrix S*{I + 274)5*. It follows 
from j5j that the classical (multi-valued) Gaussian triplet 

(5(ei), . . . , 5(e,)), (E(efc+i), . . . , B{e2k)), {B{e2k+i), • • • , ^(63^)) 

is Markovian if and only if 

[S*{I + 2A)S]i3 = [S*{I + 2A)S]u[S*{I + 2A)S]^^'[S*{I + 2A)S]23. 

Since 

[S*iI + 2A)S]rs = SliI + 2A)rsSs 

and 

[S*(J + 2A)5]l2[^*(/ + 2/l)5]2-2M5*(/ + 2A)5]23 = 5r(/ + 2A)i2(/ + 2A)2-2^(J + 2A)2353, 

the matrix 5* can be removed from the condition and we have the equivalent form 
(/ + 2A)i3 = (/ + 2A)i2(/ + 2A)^^{I + 2A)23 which means that (1, 3) element of 
(/ + 2A)~^ is 0. If the quasi-free state induced by A gives a Markov triplet, then A is 
the form of fHTj) due to Theorem 12.31 In particular, {A~^)i3 = and reference to [S] 
gives the following result. 

Theorem 3.1 Let u be a quasi-free state on CCR(7ii 07^2 ©"^^3) which is Markovian 
with respect to the decomposition Tii © 7-^2 © H3. Assume that Cj, 1 < j < 3k are unit 
vectors such that ([5^ and [5l\ ) hold. Then the classical (multi-valued) Gaussian triplet 

{B{ei), . . . ,B{ek)), {B{ek+i, ■ ■ ■ , B{e2k)), {B{e2k+i, ■ ■ ■ , B^csk)) 

is Markovian, moreover, 

{B{ei), . . . ,B{ek)) and {B{e2k+i, ■ ■ ■ , B^csk)) 

are independent. 

In a final remark we compare the Markov condition for the classical multivariate 
Gaussian triplet with the CCR case. The classical condition is A12A22 A23 = A13. The 
CCR condition can be formulated as ^12/(^22)^23 = ^13 with any continuous function 
/. (This implies immediately that A13 = 0.) Therefore, the CCR condition is much 
more restrictive. 
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